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Ti . Abstract 

o 
o 



For a model Id asymmetric double-well potential we calculated so-called 
survival probability (i.e. the probability for a particle initially localized in 



> 

Q"^ ' one well to remain there). We use a semiclassical (WKB) solution of the 

C^^ ■ Schroedinger equation. It is shown that behavior essentially depends on tran- 

sition probability, and on one dimensionless parameter A which is a ratio 
C^ ' of characteristic frequencies for low energy non-linear in- well oscillations and 

inter wells tunneling. For the potential describing a finite motion (double- 
well) one has always a regular behavior. For A ^ 1 there is well defined 
resonance pairs of levels and the survival probability has coherent oscillations 



^ ■ related to resonance splitting. However for A ^ 1 no oscillations at all for 



the survival probability, and there is almost an exponential decay with the 
characteristic time determined by Fermi golden rule. In this case one may 
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not restrict oneself to only resonance pair levels. The number of perturbed by 
tunneling levels grows proportionally to \/A (by other words instead of iso- 
lated pairs there appear the resonance regions containing the sets of strongly 
coupled levels). In the region of intermediate values of A one has a crossover 
between both limiting cases, namely the exponential decay with subsequent 

long period recurrent behavior. 
PACS numbers: 05.45.-a, 05.45. Gg. 
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I. INTRODUCTION 

Double level systems and models appear in various contexts in physics, chemistry and 
biology. The recurrent interest to the topic is related mainly with fairly rich and interesting 
physics of the systems, and with the experimental activity on several classes of systems which 
can be viewed as good physical realization of double level models (including fashionable 
quantum dots, see e.g. [|I|]). Among the possible types of behavior we will particularly be 
concerned with coherent oscillations and incoherent (dissipative like) tunneling. Our goal is 
to propose a simple mathematical model to illustrate crossover from coherent oscillations to 
dissipative tunneling (decay or relaxation), which are also related to incoherent transitions in 
multidimensional oscillator systems. In a certain sense this crossover reveals many features 
of chaotic behavior. It is a common wisdom now that classical chaos is defined as extreme 
complexity of the trajectories in phase space, with the trajectories being very sensitive to 
small changes in the initial conditions 0, 0]. As well evidently that the state vector (wave 
function) of a closed quantum system strictly speaking does not exhibit chaotic motion, as a 
consequence of the unitary nature of time evolution. But in fact since in quantum mechanics 
trajectories in the phase space can not be introduced due to Heisenberg uncertainty principle, 
the standard classical concept of the stability becomes ambiguous (see e.g. 0, 0, 0, [|^, 

We put forward a simple (but yet non-trivial) model of Id asymmetric double-well po- 
tential which can be used to describe under relatively weak assumptions a crossover from 
coherent oscillations (say mechanical behavior) to incoherent decay or dissipative tunneling 
(say ergodic behavior). The essential part of the model we will present is to illustrate this 
semiclassical quasi-chaotic behavior. In fact the illustration was made long ago by Fermi, 
Pasta, and Ulam |p. They performed computer studies of energy sharing and ergodicity 
for weakly coupled systems of A^ oscillators. Later on, the results of [^] were confirmed 
and refined (see e.g. [|1^, [|ll|])- But all these papers were devoted to systems with many 
degrees of freedom (A^ ^ 1 dimensional phase space) for the cases where the motion is 



nearly integrable and irregular in different energy regions. Level statistics for such kind of 
mixed systems (i.e. when behavior is regular and chaotic in different phase space regions) 



changes gradually from Poisson to Wigner type of distributions [0, |T3|, ||I4[. Thus these 
systems become non-integrable when the energy exceeds a certain critical value. Just on the 
contrary we will propose and investigate in Id a conservative system with time independent 
Hamiltonian which is evidently always integrable, and it does not generate classical chaos. 
For the sake of completeness let us note that the tunneling in the mixed (i.e. regular- 
chaotic) systems has been studied as well for two level systems when one of the levels 



interacts with a chaotic state [T^, ||16[ (see also review [|T^] and references therein). In the 
case of a resonance between the tunneling doublet and suitable chaotic states, the tunneling 
is enhanced (so-called chaos assisted tunneling) and has very strong resonance dependence 
on quantum numbers. Similar effects due to transverse vibrations take place for isolated 



Fermi resonances in tunneling systems [|18 



Our paper has the following structure. Section 11 contains basic equations necessary for 
our investigation. Section III is devoted to the calculation of so-called survival probability. 
We use the semiclassical approach |]TU| (see also [^ and references herein). The last section 
IV contains the summary. The appendix to our paper is devoted to the technical and 
methodical details of the calculations. 

II. ASYMMETRIC ID DOUBLE- WELL POTENTIAL 

The simple model studied in this paper consists of a quantum particle in one dimensional 
asymmetric double-well potential U{X) with one-parameter dependent shape. Using the 
tunneling distance ao and the characteristic frequency of the oscillations around the left 
minimum Qq, we can introduced the so-called semiclassical parameter 7 = mVtQa^/h ^ 1 
{m is a mass of a particle, and further we will set /i = 1 measuring energies in the units 
of frequency), which is assumed to be sufficiently large, i.e. the tunneling matrix element 
should be small in VLq scale. The choice of the model potential is dictated by the principle 



of minimal requirements. Our aim is to describe in the frame work of one universal model 
crossover from symmetric double-well potential to so-called decay potential, and to do it 
we need a parameter to make the right well (i?-well) deeper and wider than the left well 
(L-well). 

Using Qq and Oq to set corresponding scales, the model potential satisfying these minimal 
requirements can be written in the following dimensionless form 



V{x) = -x^{l-x) 



1 + ^^ 



(1) 



where V = U/{Qo'~f), and x = X/aQ. The dimensionless parameter b allows us to change the 
shape of the right well (i?-well), and to consider both limiting cases, namely a traditional 
symmetric double- well potential (for 6=1), and for 6 ^ oo a decay potential (or by other 
words to change the level spacings from Qq^ scale to zero ). In fact it can be shown (see 
below and the Appendix to the paper) that quahtatively all our results do not depend on 
the concrete form of the one parametric potential satisfying these requirements (only on the 
density of i?-states) . Behavior in both limiting cases are well known, and for 6 = 1 one has 
coherent quantum oscillations, typical for any two-level systems, while for b -^ oo there is 
a continuum spectrum of eigen states for x -^ -|-oo and one can find an ergodic behavior 
(incoherent decay). Our main goal in this section is to study crossover between both limits 
at variations of b. 

The general procedure for searching semiclassical solutions of the Schroedinger equation 
with the model potential (|l|) has a tricky point. The fact is that in the L-well we have a 
discrete eigenvalue spectrum (stationary states) while for the i?-well in the case 6 ^ 1 we 
have quasi-stationary states, which are characterized by wave functions \1/„(X) exponentially 
increased in the region of e ^ V{X). Both kind of states are defined on different sheets 



of complex energetic surfaces |T9[, and to treat both kind a states one should use different 
tools, namely, the standard quantization of the stationary states from the discrete part of 
the spectra |]T^, and proposed long ago by Zeldovich |^ for quasi- stationary states the flux 



probability conservation law, which leads to the Lorentzian envelope for spectral distribution 



functions. Unlike ||2T| in our case we get the Lorentzian envelope filled by (5-peaks of the 
final states. 



The procedure is described in the Appendix, and it includes three steps (see |T9|, ^1 
and we will use notations from [pO|]): 

• First one should find the action Wl in the classically allowed region (i.e.^^L between 
turning points) in the left well (L-well), and apply the semiclassical quantization. For 
the low energy states in the L-well it leads to the following relation 



jWl = TC 



n+- + Xr. 



TTSn , (2) 



where, integer numbers n numerate eigenvalues, Xn is determined by an exponentially 
small phase shift, and the last r.h.s. of (0) is in fact the definition for eigenvalues £„. 

Second, the same should be done for the right well (i?-well). The calculation is almost 
trivial in the limit 6^1 (when the potential (|l|) becomes strongly asymmetric) 

lWR = -fWP + nPe, (3) 

where the dimensionless energy e is counted from the bottom of the L-well, the action 



1^7 is 



(0) 

ID 

(0) 



and 



7w^r = Y^(&^-ir(&'+i), (4) 



62 + 1 
/?=^— ~6,for 6>1, (5) 



Note that the parameter (3 = 2VLq/ujii is proportional to the density of states in the 
i?-well {ujpi is the frequency of non-linear oscillations in i?-well at e = 0), and therefore 
knowing the magnitude (3 one can compute the density of states in the i?-well, which 
grows proportional to 6 for 6 ^ 1. It is convenient to rewrite (|^) - (^) in the same 
form as (|l]) 

6 



^Wr = 71 



nR + - + an + Px 



(6) 



where ur and a„ are integer and correspondingly fractional parts of the quantity 






(7) 



The physical meaning of a„ is the deviation from a resonance between the n-th level 
in the L-well and the nearest level in the i?-well. By the definition of a fractional part 

|a„| < 1/2/3. 

• And as the last step, again using the quantization rule, one can find the spectrum. 

It turns out (see Appendix) that the spectrum and the behavior of the system depends 
crucially on the parameter A = (3Rn, where 

2n+2 n+l/2 



Rn 



TT^/^ri! 



■exp{-2^WB) 



(8) 



is the (3 independent decay rate of the n-th metastable state of the L-well at 6 ^ oo {Wb is 
the action in the classically forbidden (between turning points) region). 
For A ^ 1 solving the quantization relation (|A^), one can easily find 



1 1 

^ 2 2/3 



4 
al H — ^A - a^ 



(9) 



This expression (^ determines the resonance pairs of the levels, so-called two-level systems. 
Besides from the same quantization rule ( [A2| ) we get analytically (i.e. for arbitrary values 
of A) eigenvalues for the i?-well in the vicinity of the resonance doublet 



1 1 

"+2 + 2^ 



(m - a„)2 + —A - (m - a„] 



m = ±1,±2,.... 



(10) 



These levels are numerated by the quantum number m. 

For A ^ 1 all displacements of the levels due to tunneling are small, and two-level 
system approximation is valid (i.e. there is well defined isolated resonance pairs of levels 
with splitting oc (i?„//5)^/^). The situation becomes completely different for A > 1. In the 



limit A ^ 1 we get almost equidistant spectrum of mixed L — R levels in the vicinity of the 
following values of x (see Appendix for the details) 



_ m + 1/2 - a^ 

X =^ X.nm i 



1 

vrA 



(11) 



/3 

The given above expressions (|TU|) - (|TT]) show that the number of the perturbed by 
tunneling levels grows proportionally to vA- In Fig. 1 we have shown the displacements of 
the levels perturbed by tunneling. These displacements are decreased very rapidly for the 
levels with quantum numbers larger than vA- The scales in this figure are given by the 
semiclassical parameter 7 which relates to the L-well and the barrier. Once the scales are 
fixed the i?-well is characterized by the eigenfrequency oc 1/6 at e = (or what is the same 
by the density of states or by the action Wr in the i?-well). 

Summarizing the results of this section, thus we have shown that instead of isolated 
two level systems taking place for A ^ 1, in the opposite limit A ^ 1 there appear the 
resonance regions containing the sets of strongly coupled levels. The resonance widths are 
determined by tunneling matrix elements {HI2 = co'Li^i?,exp(— 2714^5 )/47r^ = _R„//?). In spite 
of the fact that for any finite values of A (and h) we have only the discrete spectrum of real 
eigenstates, found above mixing oi L — R states very closely resembles the representation of 
quasi- stationary states in terms of eigenstates of a continuous spectrum. This behavior can 
be formulated by other words in terms of the so-called recurrence time, i.e. the characteristic 
time when the system is returned to the initial state. For a finite motion (i.e. for a finite 
value of b) the behavior of the system remains regular. The recurrence time (i.e. in the case 
merely coherent oscillation period) is proportional to I/H12 for A ^ 1, while for A ^ 1 this 
time scales as 1/ujr (as a long-period time scale). 

III. SURVIVAL PROBABILITY 

The tunneling dynamics can be characterized by the time evolution of the initially pre- 
pared localized state ^1^(0), and by the definition the survival probability of the state is 
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p(t) = i(vi>(o)i^(t))r (12) 

For the stationary states evidently P{t) = 1, while for quasi- stationary (decaying states), 
the survival probability reads 

P{t) = exp(-n) , (13) 

where T is the decay rate which should be found, and we use lj]^^ for the time scale. 

The simplest case is the coherent tunneling dynamics of two-level states. Let us consider 
the n — n' resonance region. The eigenfunctions of isolated R and L wells, \E'^, and ^^i. If 
one has the initial state 

^(0) = ^^, 
the survival probability can be easily calculated 



pifi = \ 



l + cos {2t\l-^ 




(14) 



Normalized wave functions in the L-well can be calculated trivially, and using standard 
semiclassical wave functions for the i?-well, we are in a position to compute the survival 
probability for a general case as a function of A. The results are shown in Fig. 2. 



For A <^ 1, P(t) oscillates with characteristic time scales proportional to ifj^^ = J (3/ Rn- 
In the region A ~ 1 these oscillations are strongly suppressed. The reason for the suppression 
of oscillations is related to interference of the states with energies in the resonance region. 
As a result of the interference the total probability for the system to return back from the 
i?-well is decreased, and low-frequency modulation of coherent tunneling is raised. The 
period of the modulation grows with (3, and in the limit A ^ 1 we get the dense spectrum of 
states in the i?-well, and almost exponential decay for P{t) with /^-independent relaxational 
time r oc R^^- In this case the survival probability (i.e. the probability to keep the system 
in its initial state) for the time interval <^ l/cuij decay almost exponentially with time, and 
the characteristic relaxation time r is determined by Fermi golden rule, i.e. r^^ oc HI2/UJR. 
This result is also conformed to van Hove statement [^ concerning quasi-chaotic behavior 
of semi-classical systems at time scales of the order of uJii/Hf2- 
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We can relate the phenoinenoin described above (i.e. almost vanishing probability for 
back-flow from the R to L well) to the Fermi golden rule for a transition probability 

Wf. = 2n\Hf,\'pf, (15) 

where Hif is the matrix element between the initial state Ei and the final state Ef, and pf 



is the density of final states. For our case {Hij = H^ = jRn/(3, and pf = (3/2) we get 
easily 

Wif = TrRn , 

which does not depend on pf. Therefore the Fermi golden rule corresponds to the limit 
when the back flow from the i?-well is totally suppressed due to the interference. 

The survival probability can be related also to spectral distribution of the initially local- 
ized in the L-well states. Indeed, by the definition of the spectral distribution S{E) of the 
initially prepared localized state is determined by the transition amplitudes in expansion 
over the eigenstates (\E'„,i?„): 

s{E) = j:\{nomn)\'6{E-E^^), (le) 

n 

and therefore 

/ + 00 
S{E)exp{-iEt)dE. (17) 

For ^^(0) = \E'f' the spectral distribution is a set of (5-peaks with Lorentzian envelope 



Crossover from the coherent oscillations to exponential decay occurs when the Lorentzian 
envelope begins to fill up by (5-peaks of the final states. Note that the width of the Lorentzian 
envelope ([T8|) does not depend on the final state density (see Appendix and also |2l|). We 
have shown the results of the calculation of the spectral distribution in Fig. 3. 
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IV. CONCLUSION 

Let us sum up the results of our paper. We investigated the behavior of a quantum 
particle in Id asymmetric double-well potential with one parameter dependent shape, which 
allows us to consider in the frame work of one universal model the crossover from the 
traditional symmetric double well potential to the decay one. We have shown that behavior 
essentially depends on transition probability, and on dimensionless parameter A which is 
a ratio of characteristic frequencies for low energy non-linear in-well oscillations and inter 
wells tunneling. For the potential describing a finite motion (double-well) strictly speaking 
one has always a regular behavior. For A ^ 1 there is well defined resonance pairs of levels 
and the survival probability has coherent oscillations related to resonance splitting. However 
for A ^ 1 there are no oscillations at all for the survival probability, and there is almost 
an exponential decay with the characteristic time determined by Fermi golden rule. In this 
case one may not restrict oneself to only resonance pair levels. The number of perturbed 
by tunneling levels grows proportionally to yA (by other words instead of isolated pairs 
there appear the resonance regions containing the sets of strongly coupled levels). In the 
region of intermediate values of A one has a crossover between both limiting cases, namely 
the exponential decay with subsequent long period recurrent behavior. 

However a number of remarks related to our results are in order. Many features often 
classified as evidences of quantum chaos in fact as we have illustrated in our model can occur 
for well defined states possessing only discrete energy levels. The deviation from two level 
system behavior, taking place for A ^ 1 has nothing to do with random or chaotic properties 
of the system. It means only that due to well known phenomenom of level repulsion the 
two level approximation is not adequate. Lorentzian envelope (see Fig. 3) we found arises 
from the interaction of a single level in L-well with a set of levels in the i?-well and not with 
appearance of level widths (imaginary self-energy contributions). 

One should distinguish between short-time and long-time behavior, and the boundary 
between them depends on the parameter A. Short-time returns (oc (3) are governed by 
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one or a small number of semiclassical paths, while long-time returns (oc -R^^) arise from 
interference between many paths. In the limit A ^ 1, exponential decay occurs for short- 
time dynamics, while the system remains regular for long-time scales, in contrast with 
chaotic models we discussed in the Introduction. Nevertheless the tunneling in the limit of 
A ^ 1 can induce vibrational relaxation for localized i?-levels. The relaxation appears due 
to tunneling recurrences, and results in redistribution of initial energy over all levels coupled 
with a single L-level. 

Main physical idea of our paper, namely that specific quasi-chaotic behavior is associated 
with the fact that one level in L- well in a certain condition (A ^ 1) is coupled to a set of 
almost dense levels in the /2-well, was discussed in the literature long ago p2[ (see also [pl[]), 
mainly qualitatively. Our achievement is that we alone seem to have propose the concrete 
and tractable analytical model to illustrate and to investigate explicitely and quantitatively 
this statement. 

In this respect our results are quite different from numerical investigations of billiard- 
type systems (see e.g. review article []T^), showing universal behavior of level spacings in 



finite chaotic systems. Our results (for the totally integrable Id model) demonstrate that 
level spacing distribution is not a specific feature of quantum systems with chaotic classical 
counterpart limit. Our finding of the equidistant regular level distribution is a result of the 
interaction of the single L-level with several (of the order of 10 for our particular choice of the 
parameters) /2-levels (which in own turn are regular ones). As well we should distinguish 



our model and dynamic tunneling ones p3|], [^. The latter assumed strong coupling of 
the tunneling system with an environment which destroys the coherence, whereas in our 
model the coherence is destroyed by the tunneling itself due to the high density of -R-states, 
breaking two level approximation. 

Note also at the very end of the paper that results presented here not only interested in 
their own right (at least in our opinion) but they might be directly tested experimentally 
since there are many molecular systems where investigated in the paper Id asymmetric 
potential is a reasonable model for the reality. And not only molecular systems, for instance 
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recently as a controllable two-level system, double quantum dots are proposed for realizing a 
single quantum bit in solid state systems. Experimentally [0] in these systems there observed 
two distinct regimes characterizing the nature of low-energy dynamics: 

(i) relaxational regime, when an excited-state electron population decays exponentially 
in time with a rate correctly given by Fermi golden rule; 

(ii) vibronic regime, when the population oscillates for some number of cycles before 
decaying. 

And what's more, at short times the averaged excited-state populations oscillates but 
has a decaying envelope. The similarity with the behavior we found in the paper is evident.^ 
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APPENDIX: 

The semi-classical wave function is represented in the well known WKB form 

^ = exp (iW) , 
The action W should satisfy to WKB equation 



^AU characteristics of our model are not specific only for Id case. For A ^ 1 one can expect 
similar behavior and for multidimensional systems. 

^ Equivalently it can be represented in the so-called instanton or minimum action tunneling path 
formalism pH] (see also pO]) in the form of ^ = exp{—jWE), which is more efficient for classically 
inaccessible parts of phase space. 
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and two turning points, which are boundaries of classically allowed regions, are situated near 
zeros of V{X) — e/7. 

For the asymmetric double-well potential (0) the Bohr - Sommerfeld [|19| quantization 
equations read 

t9{lWL)tgi^WR) = 4 exp{2^WB) , (A2) 

where Wb is the action in the classically forbidden region in between the turning points 
Xi,X2 in the left and right wells, and Wl^r are the coordinate independent actions in 
the classically allowed regions inside of the L (respectively R) well. Using the following 

expansion 

00 
tan z = y^^ 2z 



2 2/ 1'' 

z —IT \ m -\ — 

\ 2 , 



-1 



m=0 

one gets the almost equidistant spectrum of the mixed L — R levels, and in this condition the 



solution of (|A2|) leads to the expressions (|), ( |10|) presented in the main text of the paper. 



The time evolution of any initially prepared state can be described by a superposition of 
the eigenf unctions of the discrete and continuous spectra with time dependent phases. For 
the potential (|l|) with 6^1 the initial finite motion, i.e. the initial density distribution 

p{t) = f ' \^{X,t)\^dX (A3) 

JXi 

concentrated in the L-well at t = decreases exponentially with time 

p(t) =p(0)exp(-77t) . (A4) 

Eq. ( |A4|) signifies that the wave functions of quasi- stationary states have the form 

^„(X, t) = ^„(X) exp ii-iEn - Vn/2) t) , (A5) 

and the eigenvalues are complex and lies on the lower half-space of (e, r/) plane. The quan- 
tization of the stationary states of a discrete spectrum is performed by the requirement 



|^(X,t)p ^0, at|X| ^00. 
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This condition is impossible to impose to quasi-stationary states, since the wave functions 
\[^„(X) exponentially is increased in the region of e ^ V{X). The physically meaningful 
boundary condition as was noted first by Zeldovich |^ for quasi-stationary states can be 
written as a conservation law for the fiux probability from the L-well through the barrier. 
The difference between stationary and quasi-stationary states disappears as it should be at 

The expansion of the initially quasi-stationary state is dominated by the continuum 
spectrum eigenf unctions with the energies close to the real parts of the eigenvalues £„. 
These eigenfunctions have the form 



^k{X) 



( A{k)ct>l{X), X<Xm ^ 



(A6) 



y^lsin{kX + 6{k)), X>XmJ 

where Xm is the left turning point of the i?-well, the localized wave function 0^ is normalized 
to unity, and the phase is given 

k 

6(k) = So — arctan — , (A7) 

k — ki 

and 6q is fc-independent component, ki = y/2me^, k2 = kirjn/'^Sn- For the eigenfunctions 
with the energies e and e' close to £„ we get 

From (|A|), (^, and (^ in the limit £ - e' ^ we get 



■^^«-!\/|^4rf^- (A9) 



Expressions ( |A7| ), ( |A9| ) are valid for a continuous spectrum, for discrete levels the phase 



shift as well is governed by the probability flux from the i?-well into classically forbidden 
region, and instead of ( |A7| ) it leads 



5 = arctan \J Rnl3 , (AlO) 



and instead of (|A9|) one can easily flnd 
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Note that (|A11|) has almost the same form as ( |X9|) , although it depends on discrete energy- 



levels, and besides it has a different coefficient due to different normalization condition. 

The relation ( |A9|) shows that the probability density of the continuous spectrum eigen- 
states exhibits the Lorentzian distribution around the real part of the quasi-stationary eigen- 
values En- Expressions (|A9|) -( |Af 1| ) are equivalent to the spectral distribution (p!8|) presented 



in the main body of the paper. 

Few words concerning numerical results presented in the main text in the figures 1-3. 
The calculations have been performed to check: 

(i) semiclassical approximation for the model potential ([l|); 

(ii) the spectral distribution (|l^). 

We used the numerical diagonalization of the Hamiltonian matrix in the basis set of trial 
functions, which includes: so-called instanton wave functions of the L-well (see [^), and 
the WKB functions of i?-well. This basis was orthonormalized by using standard Schmidt 



method [^. For the L-well highly excited states near the barrier top have been also included. 
In all numerical calculations we set the value of ao (so-called defect of a resonance) as zero. 
All results presented on the figures do not depend on this particular choice. 

The numerical results confirm that Eq. (|T8|) is quiet accurate in the whole range of A 
where the transition from coherent oscillations to exponential decay occurs. Note that since 
i?-levels with the negative energy are not mixed with L-levels, and besides the resonance 
region is sufficiently narrow (i?„ = 0.01), we need not diagonalize huge matrices. For 
our purposes the diagonalization of the matrix 3000 x 3000 is more than sufficient to find 
eigenvalues in the resonance region around the n = L-level. 
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Figure Captions. 
Fig. 1 

The eigenvalues as functions of A for the zero-point level (n = 0) of the L-well. Dashed 
lines indicate the limits of A ^ 1, and A^l;7 = 10,ao = 0. 
Fig. 2 
The survival probability for different values of A and 7 = 10 : 

(a) A = 0.02 ,6 = 5 (solid hue) ; A = 0.5 , 6 = 116 (dashed line) ; 

(b) A = 0.5 , 6 = 116 (solid hue) ; A = 4.0 , 6 = 929 (dashed line) ; 

(c) A = 4.0 , 6 = 929 (sohd line) ; A = 16.0 , b = 3715 (dashed hue) . 
Fig. 3 

The spectral distribution for different values of A and 7 = 10: 

(a) A = 0.02 ,6 = 5; 

(b) A = 4.0 , 6 = 929 ; 

(c) A = 20.0 , b = 4644 . 
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